Impurity clusters and localization of nodal quasiparticles in d-wave superconductors 
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The long disputed issue of the limiting value of quasiparticle density of states p(0) = p(e — * 0) 
in a d-wave superconductor with impurities (vs its linear vanishing, po(s) oc |e|/A, near the nodal 
point e = in a pure system with the gap parameter A) is discussed. Using the technique of 
group expansions of Green functions in complexes of interacting impurities, it is shown that finite 
p(0) value is possible if the (finite) impurity perturbation V is spin-dependent (magnetic). The 
found value has a power law dependence on the impurity concentration c: p(0) oc Patc™, where pjv 
is the normal metal density of states and n > 2 is the least number of impurities in a complex 
that can localize nodal quasiparticle. This result essentially differs from the known predictions 
of self-consistent approximation: p(0) oc piv\/c/pjvA (for the unitary limit V — > oo) or p(0) oc 
(A/cV 2 ) exp(— A/cV 2 pn) (for the Born limit |V|pjv <§C 1)- We predict also existence of a narrow 
interval (mobility gap) around the Fermi energy, where all the states are localized on proper impurity 
clusters, leading to exponential suppression of low-temperature kinetics. 

PACS numbers: 71.23.-k; 71. 23. An; 71.55.-i; 74.72.-h 



INTRODUCTION 



Key objects for the low-temperature physics of high-T c 
superconductors with planar lattice structure and rf-wave 
superconducting (SC) order are quasiparticles of low en- 
ergy, e <C A -c £f (where A is the SC gap parameter and 
£f is the Fermi energy). In a clean d-wave system, they 
have "conical" energy dispersion: £k sa h\J v\q\ + ujt/f, 
if the wave vector k is close to the nodal points kj 
(|kj| = fcp, j = 1, ... ,4, Fig. Q} on the Fermi surface. 
Here q± and qi are the radial and tangential components 
of the relative quasimomentum q = k — kj and the "gap" 
velocity va is small beside the common Fermi velocity 
vp: ua/^f = A/ep =|3< 1. This defines the linear en- 
ergy dependence of quasiparticle density of states (DOS): 
Po(e) ~ pjv|e|/A {p N is the normal metal DOS), and re- 
spective power laws in temperature for thermodynamical 
quantities, in agreement with experimental data 

But this "clean" picture can fail at the lowest tem- 
peratures, dominated by the lowest energy excitations. 
By many analogies in the condensed matter, such exci- 
tations are expected to be the most sensitive to presence 
of impurities (almost inevitable in high-T c compounds). 
In particular, impurities can produce low energy reso- 
nances in the d-wave materials [2|, |3j , and even the possi- 
bility for impurity localization of quasiparticles was dis- 
cussed 0, H| • Traditionally, such effects were treated us- 
ing certain forms of single-impurity approximation, as 
T-matrix approximation || or its generalization to the 
self-consistent T-matrix approximation (SCTMA) @, 0]- 













— k ^ / 










K 









FIG. 1: Schematic of local coordinates near nodal points in 
the Brillouin zone of a rf-wave superconductor. 



Then, within the SCTMA scope, a finite limiting value 
of DOS, p(0) = lim e ^o PfeL was predicted either in the 
Born limit, IVJojv ^ 1 & n d in the unitary limit, 
I^Ipat ^> 1 [4j, |6|, for impurity perturbation V. However 
it was shown later that, for a properly defined SCTMA 
solution, p(0) should vanish (at any finite V) The 
issue of finite or zero limit for DOS is of crucial impor- 
tance for low temperature physics of real high-T c sys- 
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terns. Thus, supposing existence of p(0) ^ 0, P.A. Lee 
arrived at the prediction that the quasiparticle electric 
a.c. conductivity <j{uj) in the static limit u> — ► should 
tend to a universal value <r(0) = (e 2 /tt 2 H(3) 0], indepen- 
dent of impurity parameters at all! Similarly, a universal 
limit was predicted for the ratio of thermal conductivity 
to temperature, (k/T) t ^ — * [nk 2 , / (Shv-pv t^d)}^ -1 + (i) 
[HI Ei (where n is the number of conducting planes 
per unit cell and d the distance between them). But no 
such behavior would be possible for the above referred 
alternative of vanishing DOS at e — > 0. This alterna- 
tive was predicted in various analytical forms: sublinear 
p(e) oc lei'' with £ < 1 [l2| (also concluded from numer- 
ical solution of Bogolyubov-de Gennes equation s I l3l| ). 
linear pie) oc |e| (for non-magnetic impurities') 1 1-411 , or 
superlinear p(e) oc e 2 (for magnetic impurities) [14j and 
p(e) oc | er | / ln 2 (A/|e|) [l5|. It should be also mentioned 
that some field theory approaches to this problem even 
resulted in diverging DOS: p(s) oc 1/ [|e| ln 2 (A/e)] [T3| 



or p(e) oc exp A^J\n(K/e) /\e\ [17|. Up to the moment, 

the experimental checks for universal behavior in the low- 
frequency electric transport ID, 19[ and low-temperature 
heat transport 2^, 21 1 do not provided a definitive sup- 
port for any of these scenarios. 

It is important to notice that either the SCTMA ap- 
proach of Ref s. [j , 11, [lit an d the field-theory models 
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of Refs. 12 



17| are only justified if the relevant 
quasiparticle states are extended (non-localized). This 
requires that a quasiparticle with wavevector k have the 
mean free path £ longer than its wavelength A = 27r/fc, 
as expressed by the known Ioffe-Regel-Mott (IRM) cri- 
terion 22]. Also the derivation of the above mentioned 
universal limits fully relies on the Kubo formula inte- 
grals over the extended states (involving a certain life- 
time t = £/v for given quasiparticle velocity v). But 
in presence of localized states, the concepts of mean free 
path and lifetime do not make sense for them, and the 
Kubo formula should be reconsidered for the correspond- 
ing range of the spectrum, e.g., as for hopping transport 
in semiconductors [231 ] . Fortunately, the analysis of DOS 
can be done independently from the controversial issue 
of universal conductivity (postponing it for more detailed 
future treatment). 

As to the limiting DOS values, it was shown [24j that 
the SCTMA approach in fact admits two solutions. One 
of them, called SCTMA-1, leads to finite p(0) ^ as 
given in Refs. @, Q, but the other solution, SCTMA- 
2, yields in vanishing p(e — > 0) — * |e|/ ln 2 (A/|e|). And 
only SCTMA-2 proves to satisfy the IRM criterion at 
e — ► 0, thus appearing as the only valid SCTMA solution 
in this limit. Applying the same check to the power law 
solutions p(e) oc \e\^, one finds them to satisfy the IRM 
criterion onl y if £ > 1, which is not actually the case 
in Refs. [l2lTl3j and [lj] (for non-magnetic impurities). 
Thus, it could be thought that the ambiguity is resolved 



in favor of the superlinearly vanishing SCTMA-2 DOS. 

However the SCTMA (or field-theoretical) analysis can 
not be considered fully comprehensive for the real quasi- 
particle spectrum in a disordered system, if the localized 
states are also admitted. A more complete description 
is possible, noting that the SCTMA self-energy is only 
the first term of the so-called group expansion (GE) [l5| , 
where this term describes all the processes of quasiparti- 
cle scattering on isolated impurities. In accordance with 
the above IRM check, such processes can not produce lo- 
calized states near zero energy. But if the SCTMA contri- 
bution to DOS vanishes in this limit, the importance can 
pass to the next terms of GE, related to scattering (and 
possible localization) of quasiparticles on random groups 
(clusters) of impurities. The essential point is that their 
contribution to DOS is mostly defined by the real parts of 
Green functions which do not need self-consistency cor- 
rections and thus remain valid for the energy range of 
localized states (where the IRM criterion does not hold). 

The known approaches to impurity cluster effects 
in disordered d-wave superconductors, either numerical 
[2BI l2r3 ] and analytical [lj|, were contradictory about 
DOS and did not conclude definitely on localization. A 
practical analysis of such effects within the GE frame- 
work for an s-wave superconductor was proposed recently 
[9(, using a special algebraic technique in the limit e — > A. 
The present work develops a similar technique for a dis- 
ordered d-wave system in the limit e — > 0. We shall con- 
sider two main types of impurities, non-magnetic (NM) 
and magnetic (M), and show that only M- impurities can 
provide a finite (but different from the SCTMA-1 value) 
DOS, replacing the SCTMA-2 solution in a narrow vicin- 
ity of the (shifted) nodal energy and manifesting onset 
of localization there. The latter should produce, instead 
of universal conductivity, its exponential suppression at 
sufficiently low temperatures. 



FORMULATION OF PROBLEM 



In more detail, we use the low energy Hamiltonian 



H = ^y^k («FglT5 + «Ag2n) V'k 



1 y 



i(k-k')-p,;,t 



(i) 



p,k,k' 



where k are restricted to vicinities of nodal points kj , the 



'k>' 



k ) are composed of normal 



Nambu spinors ip^ — 
metal Fermi operators, N is the number of lattice sites 
and fi are Pauli matrices. The second term in Eq. [T] 
corresponds to the so called Lifshitz impurity model. It 
describes quasiparticle scattering on random sites p with 
concentration c = J2 P V-^ *C 1 in lattice of N sites, and 
the perturbation matrix is: 
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V = V 



T3, 
1. 



NM, 
M, 



(2) 



for specific types of impurities. We calculate the Fourier 
transform of retarded Green function (GF) matrix 



{Ml}) 



dt (3) 



and in absence of impurities it is po(s) 



_1 Imo(e), 



where g{e) = A^Tr^Gk w {2p N e/A) ln(2zA/e) [15|, 
vanishing linearly with e — ► 0. In presence of impuri- 
ties, Eqs. [3] and [5] generally relate the limit p(e — > 0) 
with the imaginary and traceful part of the self-energy 
matrix: 7 = ±TrImE k (0). If the main contribution 
to 7 comes from the GE pair term, where the share of 
A n cos k • n is negligible beside that of A\, 7 can be con- 
sidered momentum- independent. In this approximation, 
supposing also 7 <C A, we obtain 



where (...) is the quantum-statistical average with the 
Hamiltonian, Eq. [TJ and {a, b} — ab + ba is the anticom- 
mutator. The most important momentum-diagonal GF 
is presented as 



= (G k 1 - E k )- 



(4) 



where G k (e) = [e + H (v F qiT 3 + wa<?2Ti)] / (e 2 - £ k ) is 
the GF matrix of a clean system and the self-energy ma- 
trix E k is expanded into the GE series Q: 



E k = cT 



1 + c (^n cos k • n + A 2 n ) 

n^O 

1 

x(i-Ai 



(5) 



Here the T-matrix T(e) = V (1 - GVj , with local GF 

matrix G(e) = ^'^Gki describes all the scattering 
processes on a single impurity center. It can have reso- 
nance behavior near a certain energy e res < A, such that 



Re det 



l-G(e res )V 



0, 



it 



J , however at e 

simply tends to a certain real matrix T(0), to be specified 
below. 

The next to unity term in the brackets of Eq. [5] de- 
scribes scattering on pairs of impurities, separated by all 
possible lattice vectors n, while the dropped terms are 
for impurity triples and so on. The building block of 
all GE terms is the matrix i n (e) = TV" 1 J2k e ik ' n G k T 
which represents the effective (energy dependent) inter- 
action between two impurities at given separation n. Its 
zero energy limit is 



An(0) = (. 9n f 3 +/ n f 1 )T(0), 



(6) 



with real functions g n = — (Hvf/N) ^ k e lk ' n qi/£k ancl 
,/n = -(hv A /N) Ek eik ' n <?2/£k- Then the D0S reads: 



k 



(J) 



27 A 
p(Q) w PN—r In — < pN, 

7rA 7 



(8) 



and the following task is reduced to proper calculation of 
7 in function of the impurity perturbation parameters. 

Since the above mentioned matrices T(0) and A n are 
real, the imaginary part of the GE pair term is gener- 
ated by the poles of ^1 — A^j . If there is no such 

poles, one has to search for contributions to 7 from next 
order GE terms. In this course, for an Zth order GE term, 
the imaginary part is related to the poles of the inverse 
of a certain Zth degree polynomial in A n , . . . ,A ni{l _ 1)/2 
(where ni,... ,tii(i-i)/2 are all possible separations be- 
tween Z impurities). 

Generally, in the energy spectrum of a crystal with im- 
purities, one can distinguish certain intervals where DOS 
is dominated by contributions from band-like states, sin- 
gle impurities, impurity pairs, triples, etc. [271 ] . Then, for 
instance, in the pair-dominated energy interval, each dis- 
crete peak (by an impurity pair at given separation n in 
the lattice) experiences small shifts, due to the effects by 
neighbor impurities of such a pair, different in different 
parts of the system. These shifts produce a broadening 
of pair peaks, and if it is wider than the distance between 
the peaks, the resulting continuous pair-dominated spec- 
trum can be effectively described, passing from summa- 
tion in discrete n ^ in Eq. [5] to integration in con- 
tinuous r (for r > ro where r$ ~ a). Such possibility 
was shown long ago for normal electron spectrum [2a |. 
and it is even more pronounced for the superconducting 
system where the pair contribution to DOS at a given 
energy can come from multiple pair configurations (see 
Figs. 121 II below) and this multiplicity turns yet much 
greater at involving neighbor impurities into each con- 
figuration. An approach to analogous problem in s-wave 
superconductors was proposed recently Q , using an alge- 
braic isomorphism of the matrices of interaction between 
impurities in that system to common complex numbers. 



NON-MAGNETIC IMPURITIES 

Remarkably, the same isomorphism is also found for 
matrices of zero energy interaction A n between NM im- 
purities in the cZ-wave system. In this case, we have 
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|flr*| = V~ x and / r * = 0. Direct calculation of g r and f r 
shows that all possible r* lie on nodal directions, forming 
identical finite series along them (like those in Fig. [2]). 
Also we note that poles only exist at high enough band 
filling, pn£f > (l + e ) -1 an d strong enough perturbation 
parameter, V > (pn — Qas)^ 1 , as chosen in Fig. O If 
so, it is suitable to pass in the vicinity of each r* from 
integration in the components r\ , r 2 of vector r to that 
in the components s, t of the "complex" denominator: 




FIG. 2: The poles of the integrand in Eq. [9] in function of 
pair separation vector r = (x,y), are located along the nodal 
the direct space), where one of the pole conditions, f r — 0, 
holds identically. Another pole condition, \g r \ — 1/\V\ = 
\g as — 1/V|, is reached at discrete points (here at the choice 
of parameters /3 — 0.05, Vg as = 5). The plot below shows the 
related behavior of g r — 1/V for r along the nodal direction. 



explicitly: T(0) = Vr 3 , where V = V/(l - Vg as ) 
and g as ~ pn In \\/p^e-p — 1 1 is the factor of particle- 
hole asymmetry (away from half-filling: pn^f 1/2). 
Then the interaction matrices, Eq. [6l are presented 
as A n — Vg n + iVf a T 2 , thus pertaining to the gen- 
eral two-parametric family C(x,y) = x + iyr 2 with real 
x, y. This family forms an algebra with the product 
C(x, y)C(x', y') = C(xx' — yy' , yx' + xy'), isomorphic to 
that considered in Ref. [9j and to the algebra C of com- 
mon complex numbers: (x + iy)(x' + iy') — xx' — yy' + 
i(yx' + xy 1 ). By this isomorphism, the real matrix A n is 
related to a "complex number" A n = Vg n + iVf n , where 
the "imaginary unity" i corresponds to the real matrix 
i = iT2- Using such "complex" representation and the 
above mentioned passage from summation in n to inte- 
gration in r, we can write the pair contribution to 7 in 
the form: 



c 2 Im 



r>r 



dr ^ p(r) + iq(r) 
a 2 s(r)+it(r) 



(9) 



with s = 1 — V 2 (g 2 — f 2 ), t = —2V 2 g r f r and some con- 
tinuous functions p(r),q(r). Here the symbol 3? means 
the "real" (traceful) part of the "complex" integrand, 
while the common imaginary part Im is produced by its 
poles. These are attained at such separations r = r* of 
an impurity pair that s(r*) = t(r*) = 0, which requires 



Im 5Z / dsdt Jr-(s,t) 



p(s, t)s + q(s, t)t 



(10) 



and for any r* the transformation Jacobian 



Jr*{s,t) 



d(ri,r 2 ) 



d(s,t) 



is real and non-singular. Then the singularity in the de- 
nominator of Eq. [TO] is canceled by the vanishing residue 
in the numerator and, even at formal existence of poles 
in Eq. [9l they give no contribution to the zero energy 
DOS. Mathematically, this simply follows from an extra 
dimension at 2D integration, giving zero weight to the 
isolated poles. 

The above conclusion can be immediately generalized 
for any Zth order GE term (I > 3), where the integrand 
is again presented as (p + iq)/(s + it) and p,q,s,t are 
now continuous functions of Ni = 2(1 — 1) independent 
variables (components of the vectors ri, . . . ,rj_i) in the 
configurational space Si. This integrand can have sim- 
ple poles on some (Ni — 2)-dimensional surface Ai in 
Si (under easier conditions than for 1 = 2). Then the 
A'j-fold integration can be done over certain coordinates 
Hi, ... , Ujv,_2 in Ai and over the components s, t of the 
"complex" denominator in the normal plane to Af. 



Im / dui . . . du N ,^ 2 J (ui, 



, u Nl - 2 ) 



dsdt 



p(...,s,t)s + q{...,s,t)t 
s 2 + t 2 



(11) 



with a non-singular Jacobian J . The latter integral has 
no imaginary part by the same reasons as for Eq. 1101 
Thus, it can be concluded that perturbation by NM- 
impurities in a d-wave system can not produce localized 
quasiparticles of zero energy, and this directly follows 
from the indicated isomorphism of the interaction matri- 
ces to the algebra C of complex numbers. 

Moreover, the same conclusion is also valid for yet an- 
other type of NM-perturbation, due to locally perturbed 
SC order by the matrix V = Vt%. In this case, the inter- 
action matrix: 



A r 



V 



1 - V 2 g 2 a; 



[fn + Vg as g n + i(g n ~ Vg as f n ) t 2 ] , 




FIG. 3: Circular trajectories in the space of variables g T ,fr, 
corresponding to the poles of GE denominator, Eq. 1141 for M 
impurities at the choice of V = 1/pN and g as = 0.84pjv. 

pertains to the same family C{x,y) as in the above case, 
hence leading to the same absence of contribution to the 
zero energy DOS. 



MAGNETIC IMPURITIES 

However, an essential cluster contribution to zero en- 
ergy DOS can be produced by M impurities with the 
scalar (see Eq. [3]) local perturbation V = V, the respec- 
tive T-matrix in the zero energy limit being: 



OOOOOo. 




50 100 150 

k f r 



FIG. 4: In the r-plane, the trajectories of zeroes of D r form 
continuous loops within the stripes of ~ vp / (v^\8\kp) length 
and ~ l/(|<5|fcp) width along the nodal axes. The plots below 
present the pair GE term denominator D r (at the same choice 
of parameters as in Fig. [3]), seen along the stripe axis and 
along the stripe border. 



T(0) = V(l - VG)~ 



V 



1 + Vgash 
1 - V 2 9l s ' 



(12) 



Unlike the traceless T(0) for NM-impurities, this matrix 
produces a finite shift of the nodal point itself, from e = 
to e = e = cV(l — y 2 gls) 1 - As usual, this can be 
absorbed into the Fermi level position by shifting all the 
energy arguments of considered GF's, then DOS in T- 
matrix (or SCTMA) approximation will vanish at e — > 
eo (in the same way as at e — > for NM-impurities), 
fixing the relevant limit in presence of M impurities. The 
following treatment of higher order GE terms involves 
the matrix of interaction between M-impurities: 



A r = V 



(t 3 + Vgas)g r + (n - iVgg S r 2 )f r 



(13) 



Notably, it does not fit the C algebra, and though be- 
ing harder technically, this permits to expect that M- 
impuritics can effectively contribute to the zero energy 
DOS. In fact, the straightforward calculation of the cor- 
responding GE pair term leads to the general matrix ex- 
pression 



(1 A 2yl _ Np + Nlh + lN 2 T 2 + N 3 T 3 

D T 

Here A^'s in the numerator are certain functions of g r , f r 
and the denominator D is the 4th grade polynomial: 



D 



f r + (gr 



gas) y 2 



1 



fr + (gr+gas) ~ y$ 



Zeroes of D in the space of variables g r , / r form two circu- 
lar trajectories of radius 1/V centered at ±g as , as shown 
in Fig. [3] (cf. to their location in the two outer points of 
these circles, / r = 0, ±g T — g as + l/V, in the NM case). 
It is this extension of singularities, from isolated points 
to continuous trajectories, that allows finite imaginary 
part of the 2D integral 



7 



c 2 V 



i - v^gls 



-Im 



dr 7V + Vg as N 3 
a 2 D r 



(15) 



On the other hand, it is assured by the fact that the 
traceful numerator Nq + Vg as N 3 in Eq. [TSldoes not van- 
ish on these trajectories. Quantitative analysis is simpli- 
fied in the case of Vg as ~ 1, when the small parameter 
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(5 = 1 — V 2 g 2 s (such that \S\ <C 1) defines a low energy 
resonance at e res ~ AS/[VpN ln(l/<5 2 )] [2§|. In this case 
we have simply N Q + Vg as N 3 w —2g 2 /V 2 . Direct nu- 
meric calculation of the functions <? r ,/ r shows that the 
trajectories D = 0, when presented in variables r± t 2, form 
multiple loops (seen in the inset of Fig. 2]) of total number 
~ vf/(v/\\8\), each contributing by ~ 7r 2 /[(a/cF) 2 |<5|] into 
Im J g 2 /(V 2 D)dr/a 2 in Eq. [T5l Thus we arrive at the 
estimate: 7 ~ c 2 V/ 1 <5 1 , and hence to the finite residual 
DOS: 



P(eo) 



c 2 V , A|<5| 
PAr A|^ ln ^' 



(16) 



We recall that this result is impossible in the properly 
formulated self-consistent approximation [l5j |. and it is 
also in a striking difference to the SCTMA predictions, 
p(0) ~ {A/V 2 p N ) ex p(-A/cV 2 p N ) in the Born limit H 
or p(0) ~ \J cpN I A in the unitary limit The non- 
universality of this effect is manifested by its sensitivity to 
the M-perturbation parameter V, so that p(eo) is defined 
by impurity pairs only for strong enough perturbations, 
1^1 ^ (pN +3as) _1 - For weaker V, nodal quasiparticles 
will be only localized on impurity clusters of a greater 
number n > 2 (the bigger the smaller |V|), and the par- 
ticular form of Eq. \W\ would change to p(0) ~ Pnc 11 
(with some logarithmic corrections). Anyhow, a finite 
limit of DOS at zero energy is granted by the fact that in 
real high- T c systems a certain M-type perturbation can 
result even from nominally non-magnetic centers (includ- 
ing dopants) Moreover, the above considered con- 
dition |U|<? as ~ 1 does not seem very difficult, as testi- 
fied by the observation of extremely low-energy resonance 
e re s « —1.5 meV by Zn impurities in B^S^CaC^Os+a 
[33l | . It just fits the asymmetric M-resonance from Fig. [5l 
contrasting with the symmetric NM-resonance picture. 

Then the overall impurity effect in a d-wave supercon- 
ductor can be seen as a superposition (almost indepen- 
dent) of the above described effects from NM-impurities 
with perturbation parameter Vnm and concentration cnm 
and from M-impurities with perturbation parameter V 
and concentration c (supposedly c <C cnm)- An example 
of such situation is shown in Fig. [5] 

The residual DOS from pair GE term prevails within 
a certain narrow vicinity of Eq where quasiparticle states 
are all localized on properly separated impurity pairs. 
Outside this vicinity, the states are extended and rea- 
sonably described by T-matrix (or SCTMA). The tran- 
sition from localized to extended states occurs at the 
Mott mobility edges e c < Eq and e c i > £0, where GE 
and T-matrix contributions to DOS are comparable. Us- 
ing the simplest approximation for the T-matrix term: 
pie) ~ pm\e — £q\/ A, we estimate the range of localized 
states (somewhat exaggerated in the inset of Fig. [3]) as 



B 
3 



B 

en 

o 

in 
B 

o 
Q 




-0,15 -0,10 -0,05 0,00 0,05 0,10 0,15 



Energy (in A units) 



FIG. 5: Low-energy d-wave DOS p(e) in simultaneous pres- 
ence of NM-impurities (with cnm = 3% and Vnm = 1/pjv), 
producing two symmetric broad resonances, and M-impurities 
(with c = 0.03%, V — Vnm, and Vg as = 0.9), producing sin- 
gle sharp resonance at extremely low energy e res . Inset shows 
the mobility edges e c and e c / around the shifted nodal point 
eo, they separate localized states (shadowed area) with almost 
constant DOS, p(e) ~ p(eo), from band-like states whose DOS 
is close to the T-matrix value (solid line). 



£o ~ ^0 



c 2 V 



In 



A\5[ 

c 2 V 



(17) 



provided it is much smaller than the distance to the M- 
resonance: 6 C <C |eo — E res \. The same estimates for the 
mobility edges follow from the IRM breakdown condition: 
e — RcS(e) ~ ImE(e), at e » eq. 

The tendency to localization of quasiparticles can be 
generally opposed by the effects of repulsive Coulomb in- 
teraction between them and this issue was also dis- 
cussed for disordered d-wave superconductors ljj, |31|. 



These field theory treatments showed that localization 
can survive at low enough temperature. The full ac- 
count of Coulomb interactions in the present GE ap- 
proach is rather complicated technically, but a simple 
estimate follows from the overall number of (supposedly) 
localized particles within the energy range, Eq. I17[ which 



is as small as ni 



l V/A8 2 )ln 2 \AS/c 2 V\ 



Since the average distance between them ~ a/V n ' oc is 
much longer than the distance between charge carriers 
~ a/^/n ~ a/ytpNEF, the effects of Coulomb interaction 
are hopefully screened out, at least for the systems far 
enough from half-filling [32]. 

Notably, localization turns to be yet possible near the 
resonance energy, e ss E res , but this requires that the 
concentration of M-impurities surpasses a certain char- 
acteristic value c res ~ (e res /A) 2 ln(A/|e res |). In par- 
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ticular, for the choice of parameters in Fig. [5l we find 
c res ~ 3 • 10~ 4 , so that this system should be close to the 
onset of localization also in this spectrum range, where 
each localized state is associated with a single impurity 
center. 

Generally, presence of localized states near the low- 
est excitation energies in the spectrum must influence 
significantly the kinetic properties of a crystal with im- 
purities, such as electric and (electronic part of) heat 
conductivity at lowest temperatures. Taking in mind 
the above referred modification of Kubo formula for the 
energy range \e — eq\ < 5 C , their temperature depen- 
dencies, instead of reaching the universal values uq and 
[k/T](0), should rather tend to the exponential van- 
ishing: ~ exp (— Sc/ksT), at low enough temperatures: 
T <C Sc/k-g,. The latter value, at the same choice of im- 
purity perturbation parameters and typical gap A ~ 30 
meV, is estimated as ~ 0.2 K. In this context, the in- 
triguing sharp downturn of k/T, recently observed at 
temperatures < 0.3 K [11] and attributed to the low- 
temperature decoupling of phonon heat channel [35j |. 
can be otherwise considered as a possible experimental 
manifestation of the quasiparticle localization by impu- 
rity clusters. A more detailed analysis of possible non- 
universal behavior of transport properties of disordered 
rf-wave superconductors will be necessary to confirm this 
conjecture. 

CONCLUSIONS 

The quasiparticle states are considered in a d-wave su- 
perconductor with impurities beyond the (self-consistent) 
T-matrix approximation, using the techniques of group 
expansions of Green functions in complexes of interact- 
ing impurities. It is shown that, if the impurity pertur- 
bation of magnetic type is present, the indirect interac- 
tion between impurities can essentially change the quasi- 
particle spectrum near nodal points, producing strongly 
localized states of non-universal character (depending 
on the perturbation strength). Experimental check for 
possible non-universal effects in low temperature trans- 
port properties can be done, e.g., in the Zn doped 
I^S^CaC^Os+a system. 
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